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ABSTRACT 

New X-ray observatories ( Chandra and XMM-Newton) are providing a wealth of high-resolution X-ray 
spectra in which hydrogen- and heUum-like ions are usually strong features. We present results from 
a new collisional-radiative plasma code, the Astrophysical Plasma Emission Code (apec), which uses 
atomic data in the companion Astrophysical Plasma Emission Database (aped) to calculate spectral 
models for hot plasmas, aped contains the requisite atomic data such as coUisional and radiative 
rates, recombination cross sections, dielectronic recombination rates, and satellite line wavelengths. We 
compare the APEC results to other plasma codes for hydrogen- and helium-like diagnostics, and test 
the sensitivity of our results to the number of levels included in the models. We find that dielectronic 
recombination with hydrogen-like ions into high (ti = 6 — 10) principal quantum numbers affects some 
helium-like line ratios from low-lying (n — 2) transitions. 

Subject headings: atomic data — atomic processes — plasmas — radiation mechanisms: thermal — 
Xrays: general 



1. INTRODUCTION 

Modeling emission from an optically-thin collisionally- 
ionizcd hot plasma has been an on-going problem in as- 
trophysics (Cox & Tucker 1969; Cox & Daltabuit 1971; 
Mewe 1972; Landini & Monsignori Fossi 1972; Raymond 
& Smith 1977; Brickhouse, Raymond & Smith 1995). Over 
time, the codes used in these papers have been updated or 
completely rewritten; current versions are the SPEX code 
(Kaastra, Mewe & Nieuwenhuijzen 1996) and the CHIANTI 
code (Dere et al. 1997; Landi et al. 1999). The fre- 
quently used MEKAL (Mewe-Kaastra-Liedahl) code (Kaas- 
tra 1992, Liedahl, Osterheld, & Goldstein 1995) embedded 
in XSPEC (Arnaud 1996) uses data that are now in the 
SPEX code; its results are similar to SPEX. We describe 
here a new plasma emission code APEC (Atomic Plasma 
Emission Code) along with the Atomic Plasma Emission 
Database (aped). 

Our primary goal is to create plasma emission models 
that can be used to analyze data from the high resolution 
X-ray spectrometers on the Chandra and XMM-Newton 
telescopes. Due to limitations on computer speed, mem- 
ory size, and the available atomic data, early plasma codes 
included only the strongest emission lines from each ion, 
and "bundled" nearby lines, reporting only their summed 
emission at a single wavelength. In addition, the data were 
stored in the code itself, and could not be easily separated 
and studied. In contrast, we have strictly separated the 
atomic data from the code, storing the atomic data in 
FITS format files, which collectively form aped. We have 
endeavored to maintain the data in aped as close as pos- 
sible to the original form, for simplicity and ease in error 



checking. In addition, aped contains error estimates for 
many wavelengths; errors on other values in the database 
are in progress. 

Currently, aped contains data for well over a million 
lines, although many of these are too weak to be ob- 
served individually. An APEC calculation done at low elec- 
tron density (1 cm"'^), with cosmic abundances (Anders 
& Grevesse 1989) contains over 32,000 unique lines whose 
peak emissivities at temperatures between 10^ — 10^ K ex- 
ceed 10^^° photons cm'^s"^. Although the literature of 
theoretical calculations for observable atomic transitions 
is far from complete, multiple calculations for some impor- 
tant rates do exist, aped contains all the different datasets 
we have collected {e.g. the entire CHIANTI v2.0 database 
(Landi et al. 1999), as well as the data referred to in this 
paper) to allow easy comparison of different rates, and of 
the effect of different assumptions on the emissivities of 
selected lines. In a subsequent paper, we plan to present a 
complete overview of the data in aped and the emissivity 
tables over the X-ray range. Here we discuss what effect 
changing the underlying atomic data has on hydrogen-like 
(H-like) and helium- like (He-like) oxygen ions. 

2. METHOD 

APEC calculates both line and continuum emissivities for 
a hot, optically-thin plasma which is in collisional ioniza- 
tion equilibrium. Although APEC can calculate the ion- 
ization balance directly (and thus handle non-equilibrium 
conditions), aped does not yet contain all the necessary 
ionization/recombination rates, so we use tabulated values 
for the ionization balance in thermal collisional equilib- 
rium. We primarily use the ionization balance calculated 
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by Mazzotta ct al. (1998; hereafter MM98) because it is 
a recent compilation that self-consistently treats all the 
astrophysically relevant ions. APED itself contains recom- 
bination (radiative and dielectronic) rate coefficients for 
oxygen. We combined these with the ionization rates from 
MM98) to create our own self-consistent ionization balance 
for these ions and found some small differences (< 25%; 
see Figure 1(b)). We do plan to include in apec/aped 
a self-consistent calculation of the ionization balance and 
level population in the future. 

Since the term "emissivity" has a number of definitions, 
we give our definition explicitly. Similar to the terminol- 
ogy of Raymond & Smith (1977) (but in photon instead 
of energy units), the emissivity of a spectral line is the 
total number of radiative transitions per unit volume, di- 
vided by the product of the electron density Ug and the 
hydrogen (neutrals and protons) density Uh in the astro- 
physical plasma. The line emissivity therefore has units of 
photons cm'^s^^. Since the number of photons emitted is 
actually proportional to the density of the ions involved, 
this definition implicitly requires an elemental abundance 
and ionization balance for the relevant ion to be specified. 

For a given electron temperature T and electron den- 
sity rie, the level populations for each ion are calculated 
from the collisional (de-)excitation rate coefficients, the ra- 
diative transition rates, and the radiative and dielectronic 
recombination rate coefficients. Excitation-autoionization 
processes have not yet been included in APED, but these 
are not significant for the H- and He-like ions in the equi- 
librium plasmas discussed here. For the H- and He-like 
isosequcnccs, we have data for all singly-excited levels up 
to principal quantum number n = 5 with the exception of 
O VII, where we have data up to n = 10. 

For the He-like ion O VII, the collisional excitation rate 
coefficients are from Sampson, Goett & Clark (1983), Kato 
& Nakazaki (1989), and Zhang & Sampson (1987) for the 
levels up to n = 5, and from HULLAC (Bar-Shalom, 
Klapisch, & Oreg 1988; Klapisch et al. 1977) calcula- 
tions for n = 6 — 10. The radiative transition rates for 
n < 5 are from TOPbase (Fernley, Taylor & Seaton 1987), 
Derevianko & Johnson (1997), Lin, Johnson, & Dalgarno 
(1977), and from the NIST database^; for n > 5, the 
rates are again from HULLAC. The wavelengths are taken 
from Drake (1988). For the H-like ion O VIII, we use the 
collisional rate coefficients from Sampson, Goett & Clark 
(1983). We have also compared with scaled values from 
Kisielius, Berrington & Norrington (1996), but find few 
significant differences in the low-temperature region where 
the data overlap. The radiative transition rates are again 
from TOPbase and Shapiro & Breit (1959) for the two- 
photon (ls2,s^So Is^^So) rate. The wavelengths are 
from Ericsson (1977). For both ions, we use the dielec- 
tronic recombination (DR) rate coefficients and satellite 
line wavelengths from Vainstein & Safronova (1978, 1980) 
for n < 3. For n = 4, 5 we used data from Safronova, 
Vasilyev & Smith (2001). Dielectronic recombination of 
O VIII into levels n = 6 — 10 of O VII are from Safronova 
(2001, private communication). Radiative recombination 
is calculated using the Milne relation and the photoion- 
ization cross-sections of Verner & Yakovlev (1995) for the 
ground states and Clark, Cowan, & Bobrowicz (1986) for 

^ http: // physics. nist .gov / cgi-bin /AtData/main_asd 



the recombination to excited states. 

3. RESULTS 

We choose these oxygen ions for our comparisons since 
they are relatively simple and have strong emission lines in 
many astrophysical plasmas. We will concentrate on two 
line ratios: for H-like ions, Ly/?(16.020A)/Lya(18.987A), 
and for He- like ions, the frequently used G = {F + 
I)/R ratio, where F (22.098A) is the forbidden transition 
ls2s3Si ^ Is^ iSo, / (21.804A, 21. 801 A) is the sum of the 
two intercombination transitions l,s2p ■^Pi.2 Is^ ^Sq, 
and R (21.602A) is the resonance transition ls2p ^Pi 
Is^ ^Sq. Figure 1(a) shows the H-hke Ly/3/Lya line ratio 
from APEC (with the MM98 ionization balance), spex 1.1, 
and Raymond-Smith. The MEKAL results track those of 
spex 1.1, and so are not independently considered. We 
also did a "pure-APEC" calculation where the ionization 
balance was found using recombination rate coefficients 
taken from aped. However, this did not affect the O VIII 
Ly/3/Lya ratio and so it is not shown. Figure 1(b) shows 
the He-like G ratio from the same calculations plus those 
of Bautista & Kallman (2000) and the pure-APEC calcu- 
lation. The difference between the apec -I-MM98 and the 
pure-APEC calculation for O VII, and the lack of a change 
in O VIII, can be understood in terms of the excitation 
mechanisms for the lines, as will be discussed below. 

Figure 1 demonstrates two important points. First, near 
the peak emissivity of the emission lines, the code predic- 
tions agree to within 25%. Second, substantial disagree- 
ments (exceeding 50%) exist even for these simple ions 
when they are not near peak emissivity. Of course, com- 
paring three theoretical models against each other does 
not provide an accurate estimate of the model error. For 
one limiting case, we can compare with analytic results. In 
Figure 1(a), the Bethe limit of 0.14 [photon units] (Burgess 
& Tully 1978) is also plotted; this is the ratio of the elec- 
tron collisional excitation rates in the high temperature 
limit, taking into account the branching ratio of Ly/? from 
the 2>p level. In the H-like case, this limit should be a 
good estimate of the Ly/3/Lya ratio since collisional exci- 
tation is the dominant method of populating excited lev- 
els. Wc calculate that cascades from collisional excitation 
to higher levels as well as radiative recombination con- 
tribute only ~ 10% of the direct rate to each line, some- 
what more for Ly/3 than Lya. Since dielectronic recombi- 
nation cannot occur onto bare ions, it is not an issue. The 
agreement between the apec result and the Bethe limit is 
thus reassuring; as expected, apec is slightly larger due to 
cascades and radiative recombination. This also explains 
why changing the ionization balance (the pure-APEC case 
in Figure 1(a)) does not affect the O VIII Ly/J/Lya ra- 
tio. Both lines are primarily populated by direct collisions 
of electrons with O VIII, so small changes in the O VIII 
abundance cancel. 

Some Chandra and XMM-Newton measurements of the 
O VIII Lyf3/hya ratio from astrophysical plasmas already 
exist. The Chandra HETG observation of the RS CVn 
star II Peg obtained a value of 0.201 ± 0.015 for this ratio 
(Huenemoerder 2000, private communication). In addi- 
tion, both Chandra HETG and XMM-Newton RGS obser- 
vations of the oxygen-rich SNR E0102 have been made. 
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Fig. 1. — (a) O VIII Ly/3/Lyo ratio [photon units] as a function of temperature for a low electron density plasma in collisional ionization 
equilibrium. At the temperatures plotted, each line emissivity is > 10~^* ph cm^s~^, or about l/500th of the peak emissivity. The dot-dashed 
horizontal line marks the Bethe ratio. The arrow at 3.16 x 10® K marks the temperature of maximum total line emission, (b) The O VII 
G ratio [photon units], as a function of electron temperature for the same plasma codes plus data from Bautista & Kallman (2000) and a 
pure- APED calculation, including the ionization balance. For the temperature range shown, the F-|-I and R line emissivities exceed 10~^*ph 
cm*s~^. The arrow at 2 X 10® K again marks the temperature of maximum total line emission. 



and although the extended nature of the source makes 
the analysis more difficult, the O VIII Ly/3/LyQ: ratio is 
measurable. The Chandra HETG obtained 0.138 ±0.016 
(Davis ct al. 2001), while the XMM-Newton RGSl and 
RGS2 instruments obtained 0.225 ± 0.02 (Rasmusscn ct 
al. 2001). The discrepancy between these results is not yet 
understood. In addition, some laboratory measurements 
with the LLNL EBIT of oxygen Ly/S/Lya have been made. 
They find a range of values from 0.16 to 0.24 for electron 
energies above 1 keV (Gendreau 2001, private communi- 
cation). 

Although the Bethe limit of 0.14 for the Ly/3/Lya ratio 
is a good estimate if collisional excitation dominates the 
O VIII excitation, other effects could change the observed 
value. Line blending, especially of Fe XVIII and Fe XIX 
lines with O Ly/3, can increase the apparent Ly/3 emission. 
Pure radiative recombination also leads to a larger value 
of Ly/3/LyQ! in H-like ions, similarly to the well-known ef- 
fects on line ratios for He-like ions in photoionized plasmas 
(Porquet & Dubau 2000; Bautista & Kallman 2000). How- 
ever, tests with APEC show that a plasma has to be nearly 
pure 0*+ for this to occur in a hot (T > 5 x 10^ K) plasma. 
Charge exchange into excited levels has been suggested by 
Rasmussen et al. (2001) as a possible way of increasing the 
ratio in E0102, but this requires mixing of neutral mate- 
rial directly into highly ionized regions; furthermore, it is 
difficult to imagine how this could be occurring in II Peg. 
Finally, resonance scattering could also cause a change in 
the observed ratio, with the size and direction of the effect 
dependent on the spatial configuration of the plasma (e.g. 
Wood & Raymond 2000). 

Of course, the Bethe limit is only applicable if the upper 
levels of both emission lines are populated by direct excita- 
tion from the ground state, and depopulated by radiative 
transitions. In the He- like O VII system, other processes 
are more important than direct excitation at high temper- 
atures for populating the ls2s ^Si level (see Figure 2(b)). 
The differences in the O VII G ratio among the codes at 
high temperature appear to be largely due to differences in 
the dielectronic recombination rates. Figure 2 shows the 
level population mechanisms for the ls2p ^Pi level (which 



forms the resonance line) and the O VII ls2s ^Si level 
(which decays to the ground state as the forbidden line). 

The total excitation rates for these levels are equivalent to 
the line emissivities, since radiative decay to the ground 
state is the primary de-excitation process at low electron 
density for both levels. In Figure 2(a), we see that direct 
excitation dominates all other mechanisms for populating 
the ls2p ^Pi level. However, Figure 2(b) shows that for the 
ls2s '^Si level, direct excitation is important only at low 
temperatures. At the peak emissivity, the main excitation 
mechanism is cascades from higher levels. These higher 
levels are primarily populated by dielectronic recombina- 
tion, which we have checked by comparing the computed 
rates for all the processes. As a result, the Bethe limit 
cannot be used as a check on the high-temperature lim- 
its of this He-like line ratio. This can also explain the 
difference between the APEC -f MM98 and pure- APEC cal- 
culations shown in Figure 1 (b) . As a result of the differ- 
ent recombination rates (~ 30% variation), at tempera- 
tures above 2 x 10® K, the pure-APEC code has a slightly 
larger O VII population and slightly smaller O VIII pop- 
ulation. Thus the resonance line R emissivity (which is 
proportional to the O VII abundance) is increased, while 
the F-l-I emission, due largely to dielectronic recombina- 
tion from O VIII, is reduced. This leads overall to a small 
decrease in the G ratio. 

For the H- and He-like isosequences of ions with sig- 
nificant cosmic abundances (C, N, O, Ne, Mg, Al, Si, S, 
Ar, Ca, Fe, and Ni) APED includes atomic calculations for 
principal quantum numbers n < 5. We have tested the 
adequacy of this model by running APEC repeatedly using 
only the n < 2, n < 3, and the n < 4 levels of O VII, and 
find that, above the temperature of peak emissivity, our 
results are not convergent. We therefore obtained data 
for n < 10 for O VII; the collisional excitation rate co- 
efficients and radiative rates are from HULLAC and the 
DR satellite line data are from Safronova (2001, private 
communication). The O VII G ratio for each model, with 
n < 2, 4, 6, 8, and 10 is shown in Figure 3. Clearly, the con- 
vergence including all n < 4 is inadequate. By including 
all levels up to n = 10 we converge to within the 10%-20% 
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Fig. 2. — (a) Excitation to the ls2p ^Pi level as a function of electron temperature, with the ionization balance from MM98 for O VII 
included. The total (effective) rate coefficient is shown, as well as the contributions from direct electronic collisions, radiative (RR) and 
dielectronic (DR) transitions, and cascades. The dominant process at all temperatures is electron collisions, (b) Same, for excitation to the 
ls2s ^Si level. However, in this case electron collisions are significant only at low temperatures. Prom the temperature of peaJc emissivity 
and beyond, cascades from higher levels dominate. 



accuracy of the underlying atomic data. In the case of the 
R = F/I ratio, which is sensitive to the electron density 
for rie > 3 X lO^cm"^, fewer levels are required to reach 
convergence at T = 10^ K, as can be seen in Figure 3(b). 
At higher temperatures, such as T = 6 x 10^ K, the spread 
among models is slightly larger, but including only levels 
with n < 5 is still adequate. 

4. CONCLUSIONS 

This is the first in a sequence of papers on the new col- 
lection of atomic data aped and collisional plasma code 
APEC. By separating the code and data we can more 
easily test the convergence of our models and compare 
atomic data from different sources. We have shown that, 
for O VIII, using newer atomic data leads to substan- 
tially different results in the high-temperature limit for 
the Ly/3/Lya ratio, which agrees with some recent EBIT 
measurements. In addition, for He-like ions, high-n di- 
electronic recombination from the H-like ion can signifi- 



cantly affect the low-n line ratio G. Careful treatment of 
recombination to individual levels and cascades is neces- 
sary for other diagnostic line ratios; e.g. neon-like Fe XVII 
(Liedahl 2000). The methods described in this paper can 
be used to test the importance of such detailed treatment. 
Furthermore, in cases where recombination and cascades 
contribute significantly to the level population, the accu- 
racy of the ionization balance, not explicitly considered in 
this paper, will contribute to the total uncertainty in the 
model line ratio. 
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Fig. 3. — (a) The O VII G ratio at low electron density as a function of temperature (using the MM98 ionization balance), for five different 
atomic models which include levels with the primary quantum number less than or equal to 2, 4, 6, 8, and 10. By including levels up to 
n=10, the result has converged at high temperatures to the accuracy of the underlying atomic data, (b) The O VII R ratio at T = 10^ K as 
a function of electron density Ue, for the same five models. 
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